A novel insertional mutation in the connexin 46 (gap junction alpha 3) gene associated with autosomal dominant congenital cataract in a Chinese family.

PURPOSE
To identify the genetic defect associated with autosomal dominant congenital cataract (ADCC) in a Chinese family, in which 11 individuals across four generations are affected with coralliform cataract.


METHODS
Exome sequencing was performed in two of the ADCC-affected family members to scan for potential genetic defects. Sanger sequencing was used to verify these defects in the whole family.


RESULTS
By combining whole exome sequencing and Sanger sequencing, the genetic defect was revealed to be a insertion of a cytosine after coding nucleotide 1,361 (1361insC) in the gap junction alpha 3 (GJA3) gene, causing a frameshift at codon 397 (p.Ala397Glyfs×71). This frameshift mutation cosegregates with the ADCC-affected pedigree members, but is absent in unaffected relatives and 100 normal individuals.


CONCLUSIONS
A 1361 insC mutation in the C-terminus of GJA3 is found to be associated with autosomal dominant congenital coralliform cataract. This finding is similar to that of a previous publication, thus providing further evidence that the GJA3 C-terminal domain is also its mutation area, and further expanding the mutation spectrum of GJA3 in association with congenital cataract.

a novel 1361 insC (CX46A397fs) mutation in GJA3 exhibiting a difference of only 17 amino acids from 1137insC (CX46S380fs). CX46fs380 did not traffic properly and was retained within the Golgi compartment [18], resulting in a loss of function in the formation of gap and hemigap junctional channels [19]. Therefore, it is necessary to predict the function of CX46A397fs in the pathogenesis of congenital cataract.

Patient identification:
A four-generation family with autosomal dominant cataract was examined at the First Affiliated Hospital of Zhengzhou University, Henan, China. Both affected and unaffected individuals were subject to detailed ophthalmic examinations, including A-scan and B-scan ultrasonography; ocular motility assessment; and visual acuity, intraocular pressure, papillary, fundus, and slit-lamp examinations. The relevant phenotypes were documented using slit-lamp photography, and no other history of disease was examined. Eight members of the ADCC family participated in the study, and seven were diagnosed as affected. EDTAanticoagulated venous blood samples were collected from all participants. Genomic DNA was extracted from peripheral blood lymphocytes through standard procedures using Flexi-Gene DNA kits (Qiagen Inc., Hilden, Germany). This study was approved by the appropriate ethical review committees and conducted in accordance with the Declaration of Helsinki principles.
Exome sequencing: The identity of the mutant gene for congenital cataract was not initially known, so two unoperated affected individuals (III-5, IV-1) and one unaffected individual (II-3) in the ADCC family ( Figure 1) were selected for exome sequencing. The methods of exome sequencing were outlined in previous literature [20,21]. Sequencing was performed with 100bp paired-end reads on a HiSeq2000 (Illumina, San Diego, CA). Reads were aligned to the human reference genome with BWA. Variants were called with GATK and annotated with Annovar. Pathogenicity predictions for variants were obtained from Annovar and the scaleinvariant feature transform (SIFT) algorithm. Sequencing coverage depth was calculated using BEDTools and genomic coordinates provided by Illumina. Exome capture was carried out using TruSeq Exome Enrichment kit (Illumina) according to the manufacturer's protocols. Briefly, genomic Figure 1. Pedigree of family with autosomal dominant coralliform cataract. The photographs of the proband (III:5) and IV:1 were taken after birth. They were diagnosed with coralliform cataract shortly after birth. This type of cataract turbidity changes little with age and was congenital and nonprogressive. These patients denied eye trauma operation history. No other abnormalities were found except for the coralliform cataract after external ocular examination. Symbols marked with asterisks represent individuals who were analyzed. DNA samples were randomly fragmented by Covaris with a base-pair peak of 200-300 bp for the resulting fragments, and adapters were then ligated to both ends of the fragments. The adaptor-ligated templates were purified using Agencourt AMPure XP beads, and fragments with insert size 300-400 bp were excised. Extracted DNA was amplified by ligation-mediated PCR, purified and hybridized to the SureSelect Biotinylated RNA Library (BAITS) for enrichment. Hybridized fragments bound to the strepavidin beads, whereas nonhybridized fragments were washed out after 24 hr. Captured ligation-mediated PCR products were subjected to an Agilent 2100 Bioanalyzer to estimate the magnitude of enrichment. Each captured library was then loaded on a HiSeq 2000 platform, and paired-end sequencing was performed with read lengths of 100 bp, which provided at least 50× average coverage depth for each sample. Raw image files were processed by Illumina base-calling Software 1.8 for base calling with default parameters.
Mutation identification by Sanger sequencing: Per the results of exome sequencing, we used polymerase chain reaction (PCR) to amplify the mutation area of GJA3 from genomic DNA (the forward primer 5'-AGC CCT CAG CGA CCA GAT T-3' and the reverse primer 5'-CAG TCC GCC AAG CTC TAC AA-3', annealing temperature: 57 °C), and the reverse primer was used as sequencing primer. The sequencing was performed using an ABI BigDye Terminator Cycle Sequencing Kit (Applied Biosystems, CA) on an ABI PRISM 3130 DNA Analyzer (Applied Biosystems), and data were analyzed using sequence analysis software, version 3.4.1 (Applied Biosystems). Genomic DNA from all affected and unaffected members of the ADCC family, and from 100 unrelated normal controls, were analyzed. Sequence data were compared with the GJA3 reference sequence (GenBank NM_174917.2) using Sequencher 4.10.1 (Gene Codes Corp.). Nucleotide numbering reflects complementary DNA (cDNA) numbering, with +1 cor-responding to the A of the ATG translation initiation codon in the reference sequence).
Bioinformatics analysis: DNAMAN software was used to predict the functional effect of the identified (GJA3 fs397) mutation on the conservation of GJA3 protein in different species. The effects of CX46S380fs and CX46A397fs mutation on the secondary structure of GJA3 were predicted using the Garnier-Osguthorpe-Robson (GOR) method, and the alignment of these two frameshift mutation sequences was predicted using the webpage Embnet.

RESULTS
Clinical findings: After reviewing the clinical examinations and hospital records of a four-generation family diagnosed with congenital cataract (Figure 2), seven affected individuals were diagnosed with congenital coralliform cataracts, exhibiting an almost sea coral or crystalline cataract appearance. The left eye of III:5 demonstrated coralliform cataract characterized with central radial lenticular opacities resembling sea coral (Figure 2A). The right eye of IV:1 showed snowballshaped clumps of crystals (coralliform or crystalline cataracts; Figure 2B).
Gap junction alpha 3 (GJA3) mutation analysis: After all candidate genes were screened, a cytosine insertion after coding nucleotide 1361(c.1361_1362insC) was found in the GJA3 gene. This novel insertion mutation causes a frameshift immediately after codon 397 ( Figure 3A), and has been confirmed to cosegregate with the ADCC phenotype in all affected family members only. These mutations were not observed in any current database, including the HapMap and 1000 Genomes projects. Thus, the mutation is unlikely to represent single-nucleotide polymorphism.
The 1361 insC (Ala397Glyfs×71) mutation induces a significant frameshift at amino acid 397 in the GJA3 C-terminus, resulting in the mistranslation of the final 40 amino acids of the reference protein and the addition of 31 amino acids to the C-terminus of the mutant protein before an in-frame translation stop codon is detected (Figure 4). Alignment of the 71 novel amino acids in the SWISS-PROT database failed to detect significant homology with other known proteins, so the homology model of the mutant GJA3 protein is difficult to generate.
Bioinformatics analysis results: Thirteen amino acids (marked by black bars) after residue 397 of human GJA3 are predicated to be 100% homologous with those of other species, suggesting the high conservation of these amino acids in different species ( Figure 3B). The CX46A397fs mutation results in an obvious mistranslation of the final 40 (including the 13 highly conserved) amino acids in the GJA3, and would thus be expected to significantly influence protein function and induce disruption of the protein's cytoplasmic C-terminus. Using the GOR method, the predicted results for secondary structure suggest that the mutant CX46A397fs differs from CX46S380fs by only 17 amino acids. Within the 17 mutated amino acids, prediction results indicate that there is one "C" (representing random coil) replaced by "E" (beta sheet) at amino acid 383, 3 "T's" (beta turns) replaced by 3 "C's" at amino acids 387-389, and 3 "T's" along with 2 "C's" replaced by 5 "H's" (alpha helices) at amino acids 392-396. There are also 7 "C's" and one "T" that remain unchanged. However, there is a "C" to "H" variation at amino acid 397 and an "E" to "T" variation at amino acid 398; these are not included in the 17 amino acid count ( Figure 5A). The CX46A397fs and CX46S380fs are parallel starting from the 397 glycine, with the 17 amino acid mutation from 380 to 396 forming the only difference between these two mutants. The alignment of the CX46A397fs and CX46S380fs mutant protein sequences shows 96.8% identity in 466 amino acid residues ( Figure 5B).

DISCUSSION
In this study, we localized the candidate gene of a fourgeneration Chinese congenital coralliform cataract pedigree in the GJA3 locus and identified a novel insertional mutation (c.1361_1362insC) of GJA3. This insertion mutation is the second reported mutation of GJA3 to have been localized to the intracellular C-terminus, indicating that the C-terminus plays an essential role in GJA3 gene function.
It is likely that reduced intercellular communication due to human GJA3 mutations trigger the formation of cataracts, which may be explained by decreased translation and/or enhanced degradation of connexin 46 [19]. Due to the high similarity of CX46A397fs to CX46S380fs in protein sequence and secondary structure, we hypothesize that CX46A397fs may also play a similar function in the pathogenesis of congenital cataract.
CX46S380fs did not traffic properly and was retained within the Golgi compartment. The last 29 of the 87 novel amino acids generated by CX46S380fs were required for impaired trafficking and loss of function [22]. These 29 amino acids also appeared in the carboxy terminus of CX46A397fs (blanketed in black in Figure 5B). Taken together, the last 29 of the novel 71 amino acids generated by CX46A397fs are necessary for the pathogenesis of congenital cataract. In this study, the 1361 insC (p. Ala397Glyfs×71) mutation results in a 210-nucleotide addition and a late-mature translation stop codon. As indicated in the Phospho Site Plus database, no modification sites are localized in the 397 amino acid site, so this mutation is not likely to influence posttranslational modifications, including phosphorylation of the GJA3 protein. However, 13 amino acids after the 397 residue of human GJA3 are predicated to be conserved in different species. This insertional mutation results in the disruption of the 13 conserved amino acids in the intracellular GJA3 carboxy terminus, which may significantly alter GJA3 function.
In conclusion, a novel insertional mutation of the GJA3 gene was identified in a four-generation Chinese family with autosomal dominant congenital coralliform cataract. The abnormal carboxy terminus generated due to a frameshift in CX46A397fs is predicted to cause loss of function, which may be pivotal to the formation of congenital cataract.